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1. INTRODUCTION 

Laser ranging systems have developed rapidly In the past few years. 
Limitations in ranging accuracy are caused by atmospheric refraction and 
scattering. In this paper the effects of atmospheric refraction will be 
discussed. 

Atmospheric refraction will increase the optical path length to an 
orbiting satellite by a few meters when the satellite is near zenith, and by 
over ten meters when the satellite is near 10“ elevation. Saastamoinen [1] 
and Marini [2] have developed correction formulas to compensate for the 
effects of atmospheric refraction on the optical pach length. These formulas 
are convenient for correcting satellite tracking data since they require only 
surface measurements of pressure, temperature, and relative humidity. 

Marini and Murray [3] investigated the accuracy of their formula by comparing 
the results to ray traces through refractivity profiles derived from radiosonde 
balloon measurements of pressure, temperature, and relative humidity. The 
difference between the range corrections obtained by ray tracing and the 
correction formula was only a few millimeters at 10“ elevation, and 
decreased to less than one millimeter near zenith. The correction formulas 
derived by Saastamoinen and Marini, and the ray traces through the refractivity 
profiles all assume a spherically symmetric refractivity profile. The errors 
introduced by this assumption are investigated in this paper by ray tracing 
through three-dimensional profiles. The results of this investigation indicate 
that the difference between ray traces through the spherical 1 symmetric and 
three-dimensional profiles is approximately three centimeters at 10“ and 
decreases to less than one half of a centimeter at 80“ . Therefore, if the 
accuracy desired in future laser ranging systems is less than a few 
centimeters, Saastamoinen and Marini's formulas must be altered to account for 
the fact that the refractivity profile is not spherically symmetric. 
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2. MARINT AND MURRAY'S FORMULA 

Marini and Murray's formula Is based upon the theory of geometric optics. 
The optical path length Is defined as the Integral of the group index of 
refraction along the ray path 


R = / n dZ . (2-1) 

e c 8 ~ 


The ray path will lie entirely in a plane if the refractivlty is spherically 
symmetric. Using this assumption, it follows that d_l_ = dr/sin 6. This fact 
simplifies equation (2-1) to 


R 

e 



( 2 - 2 ) 


where r^^, r^, and 9 are shown in Figure 1. The equation for group refractivlty 
[3] can be written as 


where 


N 

g 


f(X) 



« 80.343 f(A) I - 11.3 Y 


9650 + ^ 


7.5(T - 273.15) 
237.3 + (T - 273.15) 

X 6.11 X 10 


(2-3) 


(2-4) 


(2-5) 


X “ wavelength of radiation (microns) 
P = atmospheric pressure (mb) 

T = temperature (°K) 

Rh = relative humidity (%) . 
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The group index of refraction is given by 


n - 1 + 10"S . (2-6) 

s s 

The desired range correction is the difference between the measured path 
length, R^, and the straight path length, R. 


AR = R - R 
e 


(2-7) 


Using equations (2-2), (2-4), and (2-7) the range correction can be represented 
as the following integral: 


AR = 10 



N 



sin 6 


dr + 



dr 

sin 0 



(2-8) 


The first term can be considered the velocity error. The bracketed term is 
the difference between the geometric lengths of the ray and 'tralght line 
paths. To evaluate these Integrals Marini and Murray [3] used the perfect gas 
law, the law of partial pressures, the hydrostatic equations, and Snell's law 
for a spherically stratified medium. They obtained the following formula for 
AR 


AR 


/(A) 

f(^,H) 


A + B 


sin E + 


B/(A + B) 
sin E + .01 


where 


f((|i,H) - 1 - .0026 cos 2<j> - .00031H 

K » 1.163 - .00968 cos - .001041^ + .00001435?^ 

A » .002357P + .000141c 

u u 

B = (1.084 X 10~®) PqTqK + (4.734 x lo"®) ^ (T' A ' /K ) 


(2-9) 


( 2 - 10 ) 

( 2 - 11 ) 

( 2 - 12 ) 

(2-13) 







I- 


■'ll' 


and 

Pq ■ surface pressure at site (mb) 

Tq » surface temperature at site (*K) 

(|i > latitude of site 

H * elevation of site above sea level (km) . 

Notice that the formula only requires surface measurements. The integrand 
of equation (2-8) was expanded in inverse powers of slnG^ to obtain 
equation (2-9). The constant .01 in formula (2-9) was derived empirically 
to compensate for the neglect of higher order terms in this expansion. 

Marini and Murray investigated the accuracy of their formula by 
comparing the results to those obtained by ray tracing through spherically 
symmetric refractivity profiles. The profiles were obtained by using 
radiosonde balloon data to calculate the group refractivity using equation (2-3) . 
The refractivity was only a function of height since one radiosonde 
balloon was used in its calculation. The results of their investigation are 
summarized in Table 1. RT^^ is the correction obtained by ray tracing 
through the spherically symmetric refractivity profiles while MM is the 
correction obtained from Marini and Murray's formula [equation (2-9)]. 


TABLE 1 

RESULTS OBTAINED BY MARINI AND MURRAY FOR RUBY LASER (X = 694 nanometers) 


Elevation 

Angle 

RT, 

X 

Mean (m) 

RTj^ - MM 

Mean (cm) 

Standard 
Deviation (cm) 

00 

o 

o 

2.47 

.07 

.04 

o 

o 

3.69 

-.1 

.07 

20® 

6.91 

-.05 

.12 

10® 

13.32 

-.08 

.49 
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3. N0NSY1®1ETRICAL REFRACTIVITY PROFILE 

In our investigation a three-dimensional refractivicy profile was 
generated from radiosonde balloons which were released simultaneously from 
three separate locations around the ranging site. Radiosonde data from 
Project Haven Hop 1 [4] was used. Figure 2 gives the locations of the 
balloon release sites, which were centered around Washington, D.C. After 
the balloons were released, they were frequently blown down range as far as 
150 kilometers. An overhead view of three typical balloon trajectories is 
pictured in Figure 3. Notice that the three trajectories are very similar. 
This indicates that at least the wind conditions at the three sites are 
relatively homogeneous. 

Three tables of pressure, temperature, and relative humidity versus 
height were generated from the radiosonde data. Each table corresponded 
to a different balloon release. An example of the data from one radiosonde 
is shown in Figure 4. The data was taken from a height equivalent to the 
site’s altitude to approximately 15 kilometers. Each table's measurements 
were normalized to standard heights to simplify the ray tracing program 
(see Appendix A). This was done by Interpolating the pressure, temperature, 
relative humidity, azimuth, and elevation between two data points of the 
balloon release. The temperature, relative humidity, azimuth, and elevation 
measurements were linearly interpolated. The azimuth and elevation at a 
particular height represent the balloon's position with respect "o the 
tracking station from which the balloon was released. Using the hydrostatic 
equation, an interpolating formula for pressure can be obtained[3] 
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Figure 4. Sample of radiosonde data. Release site “ 55. Release 
Release date “ 1/20/70. 
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H2 = height of desired pressure (km) 
= height of known pressure (km) 

= virtual temperature at (®K) 
T^2 “ virtual temperature at H2 (**K) 

= known pressure (mb) 

P2 = desired pressure (mb) . 



The virtual temperature is related to T (®K) [ 3 ] by 



(3-2) 


where e is given as before by formula (2-6) and P is the pressure in mb. 

Assuming that virtual temperature vari»>s linearly with respect to height, the 
value of T ^2 can be linearly interpolated between two data point-.. Therefore 
the desired pressure value P 2 can be found by using (3-1). Using these 

interpolation formulas, tables of pressure, temperature, and relative humidity • 

can be generated to a height as high as the final radiosonde data point. Above 

this height the pressure is assumed to be a decaying exponential, while the 

remaining measurement? are assumed to be the value of the last measured data 

point. This is the same approach that Marini and Murray used to construct 

their pressure, temperature, and relative humidity tables from radiosonde 
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data [3]. The difference between Marini and Murray and our Investigation Is 
that three tables are constructed. From these tables, the refractlvlty at 
any point along the path of each balloon trajectory can be calculated using 
(2-3) . If the refractlvlty is desired at a hel^t not equal to a standard 
height, the values of pressure, temperature, and relative humidity ate inter- 
polated. The Interpolating is done between the two standard heights closest 
to the desired height. These interpolations are identical to the Interpolations 
used in creating the tables. Using these Interpolated values the refractlvlty 
can be calculated by applying formulas (2-3), (2-4), and (2-5). 

To generate a 3-dimenslonal refractlvlty profile, we will assume the 
refractlvlty varies linearly in the horizontal direction. The refractlvlty at 
any point can be obtained by interpolating between the refractivltles calculated 
at the same height along the trajectory of each of the three balloons using 
the formula 


N(r) = N + 6-N.+ (|>-sln(0)-N, (3-3) 

— r o 9 

where 0 and 9 are the co- latitude and longitude of the vector r, * (r,8,9). The 

geometry of the problem is illustrated in Figure 5. Notice in Figure 6 that 

e-r is proportional to horizontal displacement in the north-south direction 

while 9«sin(0)’r is proportional to horizontal displacement in the east-west 

direction. Since a linear variation of group refractlvlty in the horizontal 

east-west and north-south directions was assumed, equation (3-3) follows. 

The coefficients N , N„, and N, at any height, h, are calculated by 
r 9 

solving the matrix equation 
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Figure 5. Spherical coordinate system. 
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Figure 6. Differential elements In sp'iorical coordinates. 
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(3-4) 


where and are the refractivltles measured by the three radiosonde 

ballooa's a'; the given altitude, h, and ( 02 »<l' 2 )» 

CO- latitude and longitude of the balloons at the given altitude, r Is related 
to h b> 


L 


r = h + r 


0 


(3-5) 


where r^ Is tne radius of the earth at the laser site (see Figure 1) . 

N , N„, and N vary as a function of r but are constant with respect to 0 
r 0 ^ 

and 

Ihe (gradients with respect to r, 0, and (> can be calculated by 
differentl.itlng (3-3) 


9N _ 
3r " 

3N 3N^ 3N . 

~ + 0--^ + 4>-sin(e)-^ 
3r 3r 

(3-6) 


SN 

M ■ "e ♦•cos(e).N^ 

(3-7) 


I; - 

(3-8) 


The refractlviry at any position (h,6, 4>) can now be calculated. The values 
of pressure, Lemperature, relative humidity, azimuth, and elevation angle along 
the bal'' n trajectories are obtained from the three tables at a height h. 

These values are calclated by interpolating between two standard heights. 
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The azimuth and elevation angle are converted to longitude and co-latitude 

* Equations (2-3) through (2-5) are used along with the 
Interpolated values of pressure, temperature, and relative humidity to 
obtain the group refractivities at a height h along each balloon trajectory. 
The coefficients N^, Ng, and N^are calculated using equation (3-4) at a 
height h. These coefficients are then used in equation (3-3) along with 
0 and to obtain N(h,6, <|)). By this method the three-dimensional refractivity 
profile is created. 







A 


1 
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4. ANALYSIS OF THREE-DIMENSIONAL REFRACTIVITY PROFILE AND GRADIENTS 

The nonsymnetrical aspects of group refractlvity are discussed in this 
section. If group refractlvity was spherically symmetric, the refractivities 
measured along each balloon trajectory should be equal at equivalent heights. 
Figures 7 and 8 are plots of the difference between the group refractivities 
measured along each balloon trajectory versus height. Figures 7 and 8 are 
on pages 15 and 16. Figure 7 used release sites 52, 55, and 58, while 
Figure 8 used sites 53, 55, and 57. Notice Figure 8 has a larger difference 
between the refractivities than Figure 7. 

The refractivities at heights directly above each balloon release 
station were calculated using equation (3-3). Figure 9 illustrates the 
geometry of the problem. 



Balloon Release Site 

Figure 9. Geometry of refractivity calculations directly above each balloon 
release site. 


■! 



HEIGHT (km) 



HEIGHT (km) 
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The difference between the ref ractlvlt- Lee directly abow each site versus 
height is plotted in Figures 10 and 11. Figures 10 and 11 correspond to 
the same radiosonde ballocm releases used to plot Figures 7 and 8. Notice 
that Figures 7 and 10 are similar. The reason for this is that since the 
refractivities along the three balloon trajectories (N^, N 2 > and N^) differ 
by a small quantity, the coefficients N^ and N^ calculated by martlx equaticm 

(3-4) are small. Therefore equation (3-3) for refractivlty is essentially a 
function of hel^t only (spherically synmetric) . In this case, the refractivlty 
calculated along the balloon trajectory at a height h is approximately equal 
to the refractivlty calculated at the same height directly above the balloon 
release site. Figures 8 and 10 differ significantly since N^ and N^ are large 
(N^, N^, and differ by a large amount). It is interesting to note that 
the difference between a three-dimensional and spherically symnetric refractivlty 
ray trace was larger when the balloon releases in Figure 8 (and 11) were used to 
generate the refractivlty profiles. 

The 6 and 4 gradients were also calculated directly above each release 
site using (3-7) and (3-8) and plotted in Figures 12 and 13 along with the 
group refractivlty. If the refractivlty profile was spherically symmetric 
these gradients would be zero. Equation (3-3) was used to calculate the 
refractivlty. Note that both plots of group refractivlty follow an exponential 
curve. The e and ^ gradients in Figure 13 have the same general tendencies, 
while the 6 and 4i gradients in Figure 12 are generally unrelated. In both 
figures, though the 6 and 4> gradients are not zero. This Indicates that 
refractivlty is not spherically symmetric. The question then is how do 
these asymmetries affect the accuracy of Marini and Murray's range correction 


formula? 



HEIGHT (km) 



Figure 10. difference between group refractivlties calculated directly 
above each site versus height. Site 1 ■ 52. Site 2 * 55. 
Site 3 - 58. Time - 6:30. Date - 1/27/70. (N 


-- - <v-v- 




- V 






HEIGHT (km) 





REFRACTIVITY AND GRADIENTS 


Figure 12. Group refractivities and partial derivatives with respect to 

9 and <t> calculated directly above site 52. Auxiliary sites ■> 5! 
and 58. Time ■ 6:30. Date ■ 1/27/70. Partial derivatives ( 9, t) 
have a dimension of (radians)”^. 



HEIGHT (km) 




REFRACTIVITY AND GRADIENTS 


Figure 13. Group refractlvlties and partial derivatives with respect to 
0 and (|) calculated directly above site 53. Auxiliary sites = 
and 57. Time ■ 21; 30. Date • 1/19/70. Partial derivatives 
have a dimension of (radians)”^. 
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5. RESULTS FBfHf TBREE-DIMEIISIONAL REFRACTIVITT PROFILE RAY TRACING 


The error introduced Into leser ranging Beasureaents by asstming the 
refractivity profile is spherically synantric is investigated in this 
section. Ray traces through three-dinensional aiui spherically synaetrlc 
refractivity profiles were coapared to evaluate this error. Appendix 1 
gives the details of the ray tracing procedure. In creating the three- 
dinensional refractivity profile, radicMonde data froa Project Haven Hop 1 
was used. One balloon release site was assuaed to be the ranging site. 

The training tuo release sites are located approxiaately 100 kiloaeters from 
the ranging site. Figure 14 Illustrates the approximate locations of one 
set of three release sites. The three-dimensional ray trace, RT^, was 
dependent on azionith angle since the refractivity varied linearly in the 





Figure 14. Locations of one set of release sites. 
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horisontal direction. For the ezivuth angles were selected in 10° incre- 
aents frtm the beginning azlauth. to the final aziauth, The azimuth 

angles formed by the laser site and remaining two release sites, and A 2 , are 
used to calculate the beginning and final azimuth angles. is the closest 
azlBRith angle less than A^, which is an even multiple of 10°. ^2 ^s the closest 

azimuth greater than A^ which is an even multiple of 10°. See Figure 14 for 
beginning and final azimuth angles for one set of three balloon release sites. The 
spherically symmetric refractivity profile was calculated from the balloon released 
from the ranging site. RT^ was independent of azimuth since refractivity was only 
a function of height. Each spherically syasietrlc ray trace, therefore, was com- 
pared to three-dimensional ray traces at different azimuth angles. Figures 15 
through 21 represent graphs of this comparison. The data points (RT^-RT^), which 
were calculated every 10°, were connected by straight lines. In Figure 3, the 
balloon released from the ranging site was blown at an azimuth angle of approxi- 
mately 90°. When the same balloon releases as illustrated in Figure 3 were used, 
the difference betwee' RT^-RT^^ was small at an azimuth of 90° (See Figure IS for 
these results). This is expected since the 6 and ^ along the ray path are similar 
to the 6 and ^ of the balloon trajectory at a given height. 

The histograms in Figures 21, 22, 23, and 24 we’-e obtained from 428 com- 
parisons of ray traces through three dimensional and spherically symmetric 
refractivity profiles. The 428 ray traces at each elevation angle were obtained 
from 17 different sets of 3 balloon releases. For each set of releases, rays were 
traced at 8 to 9 azimuth angles from each of the three release sites. All refrac- 
tivity profiles were generated for a wavelength of 530 nanometers, which is the 
wavelength for the doubled Yag laser. The actual data for these 428 ray traces is 
in Apnendix B pages 67 to 134. It is interesting to compare the results obtained 
in this investigation to Marini and Murray's results (Table 1). At an elevation 
angle of 10° the standard deviation of RT^-RT^^ is 287 centimeters. Marini and 
Murray calculated the standard deviation of the difference between their formula and 
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Figure 15. Difference betv/een three-dimensional refractlvity ray trace 
(RT ) and spherically symmetric refractlvity ray trace (RTj^) 
versus azimuth. Laser site ■ 52. Aux. sites “ 53 and 56. 
Time - 17:30. Date - 2/13/70. 
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Pi<;ure 16. Di'^i'erence between three-dimensional refractlvity ray trace 
(RT_) and spherically syimnetric refractlvity ray trace (RT^) 
versus azimuth. Laser site ■ 55. Aux. sites ■ 53 and 57. 
Time - 21:30. Date - 1/19/70. 
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Figure 17. Difference between three-dimensional refractlvlty ray trace 
(RT-) and spherically symmetric refractivity ray trace (RT ) 
vetsus azimuth. Laser site “ 53. Aux. sites “ 52 and 56. 
Time - 17:30. Date - 2/13/70. 
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Figure 18. Difference between three-dimensional refractivity ray trace 
(RT-) and spherically symmetric refractivity ray trace (RTj^) 
versus azimuth. Laser site <■ 56. Aux. sites ■ 52 and 53. 


Time = 19:30. Date - 2/13/70. 
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Figure 19. Difference between three-dimensional refractivity ray trace 
(RT ) and spherically synmetrlc refractivity ray trace (RT^^) 
versus azimuth. Laser site ■ 56. Aux. sites ■ 52 and 53. 
Time = 17:30. Date - 2/13/70. 
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Figure 20. Difference between three-dimensional refractlvlty ray trace 
(RT^) and spherically synnetrlc refractlvlty ray trace (RT^) 
versus azimuth. Laser site » 56. Atix. sites 52 and 54. 
Time - 17:30. Date - 1/15/70. 
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Figure 21. Histogram of RT 3 -RT 1 data for E - 10® and wavelength « 530 
nanometers. Sample size « 428. 
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Figure 22. Histogram of RT 3 -RT 1 data for E ■ 20® and wavelength ■ 530 
nanometers. Sample size ■ 428. 
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Figure 23. 
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Histogram of RT 3 -RT 1 data for E 
nanometers. Sample size ^ 428. 
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Figure 24. Histogram of RT 3 -RT 1 data for E * 80® »nd wavelength * 530 
nanometers. Sample size » 428. 





ray tracing through a spherically synmetric profile to be only .49 

centiaeters at E " 10^, The standard deviations of RT^-RT^ obtained In this 
investigation are larger than Marini and Murray's standard deviations of 
RT^-Mf by a'eost a factor of 10 at elevation angles of 20°, 40°, and 80°. 

The distributions in Figures 25, 26, 27, and 28 were calculated from the 
same data used for the histograms (N > 428). The percentage of RT^-RT^ values 
greater than a value can be found from these distributions. The dashed line 
represents a theoretical Gaussian distribution with zero mean and standard 
deviation equal to the standard deviation calculated at each elevation angle. 
Notice th.~ 84Z of the RT^-RT^ values at E - 10° were greater than Marini 
and Murray's standard deviation at 10°. The other elevation angles calculated 
also had high percentages greater than Marini and Murray's standard deviation 
(85Z @ E - 20°, 85Z ? E - 40°, 91% @ E - 80°). 

The spherically synnetric reiractivlty ray traces which we evaluated were 
also compared to Marini and Murray's formula, MM. Table 2 summarizes these 
results along with the results obtained for RT^-RTj^ and RT^-MM. For our 
data the mean values of RT^-MM at different elevation angles are much larger 
than Marini and Murray's data (summarized in Table 1). We believe this can 
be partially accounted for by the fact that our radiosonde data was taken 
only during January and February of 1970 and in only one locality, the 
Washington, D. C. area. Marini and Murray used data from more than one 
locality and over longer periods of time [3]. The standard deviations between 
RT^ and MM compare very well with Marini and Murray's results. Marini and 
Murray's results were obtained for a ruby laser at a wavelength of 694 nano- 
meters. To obtain the mean between RT^ and MM the means of RT^-RT^ and RTj^-MM 
can just be added. The standard deviation of RT^-MM can be found from the 
standard deviations of RT^-RT^^ and RTj^-MM by assuming these values are 
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Figure 26. Distribution of RT 3 -RT 1 data for E - 20* and wavelength - *^30 
nanometers. Sample size ■ 428. (- - -) represents Gaussian 

distribution of zero mean and standard deviation ■ 1.16 centimeters. 
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Figure 27. Distribution of RT 3 -RTX data for E » 40® and wavelength * 530 
nanometers. Sample size * 428. (- - -) represents Gaussian 

distribution of zero mean and standard deviation « .63 centimeters 
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Figure 23. Distribution of RT 3 -RTX data for E - 80® and wavelength • 530 
nanometers. Sample size ■ 428. (- - -) represents Gaussian 

distribution of zero mean and standard deviation ■ .44 centimeters 
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TABLE 2 

SUMMARY OF RTj-RTj^, AND RT^-MM RESULTS 



RTj^-MM 

RTj-RTj^ 

RT -MM 


N-51 

N>428 


N-428 

E 

Mean 

(cm) 

Standard 

Deviation 

(cm) 

Mean 

(cm) 

Standard 

Deviation 

(cm) 

Mean 

(cm) 

Standard 

Deviation 

(cm) 

80 

.10 

.06 

-.05 

.44 

.05 

.44 

40 

.15 

.09 

1 

• 

o 

.63 

.09 

.64 

20 

.37 

00 

• 

1 

• 

o 

1.16 

o 

4 

1.17 

10 

.69 

.49 

1 

• 

o 

2.87 

00 

• 

2.91 


X =* 530 nanometers 


statistically IndepeL Using the abbreviation STD for standard deviation 

equation (5-1) results. 

STD(RT^-MM) = (STD^(RT^-RTj^) + STD^(RTj^-MM))'^ (5-1) 

Equation (5-1) was used to calculate standard deviations of RT^-MM in Table 2. 

We believe that Marini and Murray's results are valid when the refractivity 
is spherically sytnnetrlc. Saastanolnen [1] estimated that the departure of 
the atmosphere for spherical symmetry would introduce at most 1 to 2 centimeters 
error into the range correction formula at 10° elevation. However our data for 
10° elevation (Figure 26) Indicates that 68% of the values (errors introduced 
by nonsymmetrical refractivity profile) are greater than 1 centimeter and 38% 
are greater than 2 centimeters. 

The dependence of RT^-RTj^ on wavelength was examined by comparing a small 
sample (N*42) of ray traces at 530 and 353 nanometers. The wavelength for 
the tripled Yag laser is 353 nanometers. The data obtained from this wavelength 
comparison is in Appendix B, pages 135 to 144 . Notice from equation (2-3) 




that group refractivlty la inveraaly ralatad to wavelength. Figures 29, 30, 
31, 32, and 33 are plots of this coaparlson. The difference between the two 


40 


± 1 . 


frequencies is snail when RT 3 -RTj^ is small (Figure 33) and 1: when RT^-RT^ 

is large (Figure 30). Notice though, that the value of |RT^-RT^| at a lower 
wavelength is always greater than at a hl^er wavelength. This is expected 
since the refractivlty Increases with decreasing wavelength (equation 2-3). 
Figure 33 gives a good illustration of this fact. The results of frequency 
dependence are sugnarlxed in Table 3. 

TABLE 3 

WAVELENGTH (XMPARISCM OF RT^-RT^ 



RTj - RT^ 

N - 42 

A • 353 nanometers 

RT^ - RTj^ 

N - 42 

A ■ 530 nanometers 

E 

Neao (cm) 

Standard 
Deviation (cm) 

Hean (cm) 

Standard 
Deviation (cm) 

80 

-.25 

.42 

CM 

a 

.39 

40 

-.22 

.63 

.20 

.58 

20 

.02 

l.?0 

.02 

1.12 

10 

1.32 

2.55 

1.23 

2.35 , 
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Figure 30. Comparison of RT 3 -RT 1 at X - 353 and 530 nanometers 
E a 10°. Laser site ■> 56. Aux. sites “ 52 and 55. 
Time - 21:30. Date - 1/21/70. 
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Comparison of RTo-KTj^ at X = 353 and 530 nanometers 
E * 10®. Laser site *= 58. Aux. sites = 52 and 55. 
Time = 17:30. Date « 1/19/70. 
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Figure 33. Comparison of RT 3 -RT 2 at \ =353 and 530 nanometers 
E =10*. Laser site = 56. Aux. sites = 52 and 55. 
Time = 21:30. Date = 1/30/70. 
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6. ACCURACY OF RADIOSONDE BALLOON DATA 

The results in Sections 4 and 5 indicate that the asymmetries In atmos- 
pheric refractlvlty can Introduce substantial error Into the surface 
correction formulas. This Is true even though the asymaetries appear to be 
small. The variations In refractlvlty at a given height measured by the three 
radiosondes along the balloon trajectories and directly above the release 
sites are only on the order of 1 to 2 percent (Figures 7, 8, 10, and 11). 
Assured pressure differences are typically a few millibars and temperature 
differences are typically two degrees or less. These small differences 
are not substantially greater than the accuracy of the radiosonde measure- 
ments. The question then arises, "Are the effects we observed in Sections 4 
and 5 due to actual asymmetries in the atmosphere or are they caused by errors 
in the radiosonde measurements?" To answer this question we Introduced ran 'om 
errors into some selected samples of the ..ctdlosonde data and obser\’ed the 
effects on the one-dimensional and three-dimensional ray trace results. 

Typical radiosonde measurement errors [5] are summarized in Table 4. 

These errors were Introduced by adding statistically Independent random 
numbers to the radiosonde measurements. The random numbers were Gaussian 
distributed with a zero mean and standard deviation equal to the root mean 
square error of the corresponding measurement. The temperature error had 
a standard deviation of .7 (°C) . The pressure error had a standard deviation 
which was a function of height (Example from Table 4; height = 35,000 feet 
then the standard deviation ■ 1 (mb.)). The relative humidity standard 
deviation was a function of temperature. If the tenperature was below -40 °C 
a standard deviation of 10 (%) was assumed. The temperature and pressure 
errors were just added to each measurement, while the relative humidity 
error was a percentage error. This procedure of calculating the random errors 
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TABLE 4 

STANDARD ACCURACIES OF RADIOSONDE BALLOON MEASUREMENTS [5] 


Parameter' 

Operating Range 

1 

Accuracy 

(Root Mean Square) 

Temperature 

Surface - 120,000 feet 

0.7 °C 

Relative 

o 

T > 0 C 

5% 

Humidity 




0 >T > -40 C 

10% 


T < -40°C 

Questionable 

Pressure 

Surface - 10,000 feet 

. 7 mb . 


10,000 - 20,000 feet 

1.0 mb. 


20,000 - 30,000 feet 

1.2 mb. 


30,000 - 40,000 feet 

1.0 mb. 


40,000 - 50,000 feet 

.7 mb. 

1 

50,000 - 60,000 feet 



.55 mb. 


was done for each measurement. Each Gaussian error calculated was Independent 
of any other error. 

Comparisons were made between ray traces through refractivlty profiles 
generated by error free data (original radiosonde data) and error data 
(Gaussian errors Introduced). The error data consisted of two types, one 
with pressure only errors, the other with temperature, pressure, and relative 
humidity errors. Table 5 summarizes the results obtained from these ray 
traces. STD(RT 2 -RTj^) was obtained from ray traces through refractivlty 
profiles generated by error free data, while STD(RT'^-RTp was calculated 
from ray traces through refractivlty profiles generated by error data 
(pressure, temperature, and relative humidity errors). The sample size 
for these standard deviations was 80. Notice that the ratio of STD(RT^-RTj^) 
to STD(RT’^-RTp is 5 at E = 10°. This ratio can be regarded as the effective 



TABLE 5 


RESULTS OBTAINED FROM INVESTIGATION OF TYPICAL ERRORS BEING INTRODUCED INTO RADIOSONDE BALLOON DATA 



STANDARD DEVIATION (cm) 

.... . .. - - - - - 


RT^-RTj^ 

(RT^-RT^) - (;. ^-RTp 

RT^-RT^ 

RT^-RT^ 

MM-RTj^ 

MM"-RT^' 

RTj-RT^ 

MM-MM" 

E-80° 

0.26 

0.19 

0.24 

0.11 

0.08 

0.07 

0.24 

0.08 

E=40° 

0.A7 

0.23 

0.31 

0.20 

0.11 

0.12 

0.30 

0.12 

E=20° 

1.09 

0.33 

0.55 

0.27 

0.26 

0.25 

0.44 

0.23 

E=10° 

3.52 

0.72 

1.05 

0.45 

0.63 

0.71 

0.58 

0.50 

N 

80 

80 

80 

40 

12 

12 

12 

12 


RT'^ = 3-dimenslonal refractlvity ray trace with typical errors in pressure, temperature and 
relative humidity. 

RT^ = 3-dimensional refractlvity ray trace with typical errors In pressure. 

RT’^ = Spherically synmet^ic refractlvity ray trace with typical errors In pressure, temperature, 
and relative humidity. 

RTj = Spherically symmetric refractivity ray trace with typical errors in pressure. 

mm" = Marini and Murray's formula correction with typical errors in pressure, temperature, and 
relative humidity. 

= Sample Size 


N 


4S 

00 




"signal-to-noise ratio" (SNR). This "SNR,'' though, decreases to 1.38 0 E - 80°. 
The standard -’aviations of RT^-RT'^ and RT^-RT^ were calculated to examine 
what percentage of the change ..a ray traces was caused by pressure errors. 

The indications are that errors in pressure can account for approximately 50% 
of the total change in ray traces. The standard deviation of MM~RT^ is approxi- 
mately equal to STD (MM"-RT'p . This means Marini and Murray's formula is 
somewhat self-compensating when random errors are Introduced into the radiosonde 
measurements. It is interesting to compare the standard deviations of RT^^-RT^ 
and RT^-RTy At an elevation angle of 10° STDCRT^-RT'^) is greatei than 
STD(RTj^-RTp . This can be accounted for by the fact that for RT'^ random 
errors were introduced into three balloon releases which were then used to 
calculate the refractivity along the ray path, while the refractivity along 
the ray path for RT'^ was calculated us-ing only one. At 80° STDCRT^-RT'p is 
equal to STD(RTj^-RT'p because the refractivity is primarily a function of the 
refractivity measured by the balloon released from the laser site. The standard 
deviation of MM-MM" represents the magnitude of the error introduced into 
Marini and Murray's formula by surface measurement errors. It is interesting 
to note that STD (MM-MM") is approximately equal to STD(MM-RT^) at each eleva- 
tion angle. Figures 34 through 39 are graphs of RT^-RT^^, RT^-RT'^ versus 
azimuth. The "SNR" is similar between Figures 34 and 35. Notice how RT'^-RT" 
varies around RT^-RT^ in Figure 39. 

At an elevation angle of 10° the change in the ray trace results caused by 
errors in radiosonde data (STD(RT”-RT'^)j is only 20% of the total range cor»-ectinn 
error (SID(RT^-PT^)) , These results, indicate that Gaussi n errors in radiosonde 
measurements are not the primary source of the refractivity g;-adients which 
were observed In Sections 4 and 5. 
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Figure 34. Plots of RT 3 -RT] without random errors , with pressure 

errors (---), and with pressure, tempe ture, and relative 

humidity errors (-•-•-). E « 10®. Laser site ■56. g 

Aux. sites ■ 52 and 54. Time ■ 23:30. Date ■ 1/21/70. 
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Figure 35. Plots of RT3-RTj^ without random errors ( ), with pressure 

errors (---), and with pressure, temperature, and relative 
humidity errors (-•-•-). E ■ 10®. Laser site ■ 52. 

Aux. sites = 55 and 5h. Time = 21:30. Date * 1/21/70. 
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Figure 36. Plots of RT 3 -RTj^ without random errors ( ), with pressure 

errors (---), and with pressure, temperature, and relative 
humidity errors (-•-*-). E * 20*. Laser site ■ 56. 

Aux. sites ” 52 and 54. Time » 23:30. Date ■ 1/21/70. 
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Figure 37. Pious of RT 3 -RTX without random errors ( ) , with pressure 

errors (---), and with pressure, temperature, and relative 
humidity errors (-•-•-). E ■ 20*. Laser site ■ 52. 

Aux. sites - 55 and 56. Time * 23:30. Date ■ 1/21/70. 
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Figure 38. Plots of RT3-RT1 without random errors ( ), with pressure 

errors (- - -), and with pressure, temperature, and relative 
humidity errors (-•-•-). E • 40*. Laser site • 52. 

Aux. sites * 55 and 56. Time * 23:30. Date “ 1/21/70. 








:iMUTH 


; random errors ( ), with pressure 

i pressure, temperature, and relative 
E ■ 80®. Laser site » 52. 

Time - 23:30. Date « 1/21/70. 





7. SURFACE CORRECTION FORMULA 

In tills section a surface correction foranila will be developed to 
partially conqiensate for the asyanetries in the refrartlvlty. Using Figure 
40, formula (7-1) can be obtained for the refractivity at any value of h and 
a. Figure 40 is on page 57. This formula for group refractivity assumes 
refractivity varies linearly in the horizontal direction 

H,(h) - ^ H^(h)|a) (7-1) 

where ^ is the angular separation between the two sites, and Ng^ (h) and 

(h) are the group refractivities directly above their respective sites. 
Using Figure 41 and approximating 6 by E and y by 90^, the chord distance 
betwecul the two surxace points can be calculated. Assuming the chord distance 



Figure 41. Geometry for approximations made. 
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is equivalent to the arc distance, o can be approximated by h/(r^* tan E). 


Applying this approximation rearranging terms simplifies (7>1) to 


H (h) - N (h) + 

O 


\ (h) - H (h) 
*2 *1 

r ** 1 

* 

L'o !J 


(7-2) 


uhere 


“ nominal earth radius ” 6378 km 


Equation (7-2) for group refractivlty can now be substituted into integral 
(2-8) to obtain 


0'"/ 

'o 
r _ 





» (r) 

r 

*2 *1 

dr + lO"^ 

f 

tan E 

i(i sin 6 

^ J 

^ sin 9 


/ 


dr 


sin 6 


- R 


(7-3) 


Notice that the last two terms are just the integrals which Marini and Murray 
evaluated to obtain their correction formula for site 1, MM^^. The first term 
is an error term which compensates for linear variations in refractivlty in 
the horizontal direction. Thus we can write 


AR 


in 


-6 


r.v'' E 

iJ 


'1 “g 




sin 6 


dr - / -J— 


sin 6 


r dr 


+ MM, 


(7-4) 


Using Snell's Law i:i] and expanding in Inverse powers of sinCe^) gives 


sin 6 


10-^ f 


^ ~ 10 
dr - 


-6 


-6 


sin 6 


0 
-12 


/ hN„(h) dh i / h^^(h) dh 

“ sin-* Op KO * 


- / hNg(h)(NQ - N(h)] dh 


(7-5) 


where 


h - r - to . 
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N is the phase refractivlty, and is given at optical frequencies [3] by 

s-3 


f • 

N - 77.6 ^ 1+2i52_5J^ 

‘ X 

V J 


(7-6) 


An approximate evaluation of the first term of the expansion in equation (7-5) 
was calculated by Marini and Muriay [3] 


10*^ / hNg(h) dh - f(A) r^d.OSA x lO"®) PqT^K 


(7-7) 


where F^, T^, and K are as defined In an earlier section. The remaining 
terms of the expansion will be neglected since their contribution is only a 
small percentage of the final correction (millimeters compared to centimeters). 
Combining Equations (7-A) , (7-5) and (7-7) gives 


i|> tan E sin E 


M . (1.084 X 10 °) 




(7-8) 


where MM^ is the Marini and Murray correction for site 1 (equation (2-9)) 
and * function of the surface measurements at sites 1 

and 2 respectively. These measurements are taken at identical heights above 
sea level. The fact that 0^ can be approximated by E was also used in 
equation (78). Noting that where d is the arc distance between 

the two sites. Equation (7-8) becomes 


AR - MM,+ 


1 tan E sin E 


(7-9) 


where 


r 
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Foraula (7-9) was used to correct a small number (N • 10) of laser ranging 
measurements. This was done by selecting the azimuth angle for RT^ so that 
the laser beam propagated directly over one of the other release sites. The 
first two radiosonde data points of this site were used to Interpolate the 
pressure, temperature, and relative humidity at the same height as the laser 
site. In all cases the height of the laser site was greater than the remaining 
two release sites. These values were used to calculate 
(7-8). calculated from the first radiosonde data point of the 

balloon released from the laser site. The stancard deviation of the difference 
between Marini and Murray's formula and RT^ was 2.43 cm and the mean was .31 cm 
for the 10 ray traces. When equation (7-9) was used the standard deviation 
was reduced to .69 cm and the mean was .27 cm. An elevation angle of lo'^ was 
used because the error In Marini and Murray's formula Is largest at this 
angle. Figure 42 Is a graph of C/(sln£ • tanE) with the circles denoting the 
values obtained at different elevation angles by ray tracing through a three- 
dimensional refractlvlty profile. The constant C was calculated so that at 
£■>10*^ the curve C/(slnE • tanE) matched the first data point. C was not 
calculated from equation (7-9) because a second set of surface measurements 
at a point beneath the laser beam trajectory was not available. The close 
match between the theoretical curve and the ray trace data points Indicates 
that the error closely follows the l/(slnE • tanE) dependence predicted by 
the theory. 






8. CONCLUSIONS 


Numerous correction formulas have been developed to compensate laser 
ranging data for atmospheric refraction. The formulas developed by Saastamolnen 
(1) and Marini (2) are particularly convenient for use in correcting satellite 
tracking data because they use surface measurements of pressure, temperature 
and humidity to derive the range correction. The formulas were derived by 
assuming that the refractive index profile in the atmosphere is spherically 
synnetrlc. 

We have investigated the accuracy of Marini's formvla by comparing It 
with corrections obtained by ray-tracing through 3-dimensional refractivity 
profiles. The 3-dlmensional profiles were generated using measurements 
obtained from radiosonde balloons which were released simultaneously from 
three separate locations aroxmd tie ranging site. The refractivity was assumed 
to vary linearly in the horizontal direction. Also, the ray tracing method was 
based In part on the hydrostatic equation because the heights which appear In 
the radiosonde profiles were calculated from the measured pressures, temperatures 
and relative humidities using the hydrostatic equation. Our results Indicate 
that the departure of the refractivity profiles from spherical symmetry introduce 
approximately 3 centimeters error Into Marini's correction formula at 10° 
elevation. The error decreases to a few millimeters near zenith. For a given 
3-diniensional profile, the error Is a function of azimuth and elevation angles 
and can vary over a range of more than 5 centimeters for a 90° change in 
azimuth. The error decreases as the wavelength of the laser is Increased. 

A surface correction formula which compensates for horizontal variations 
in refractivity has been developed. The formula requires the pressure, 
temperature and relative humidity at the laser site and at a second location 
beneath the laser beam trajectory. Although the formula has been tested on 
only a small sample of data, the results are encouraging. 
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APPENDIX A 
RAY TRACING PROGRAM 

The first step in the ray tracing program is to read the three sets 
of radiosonde measureme;:ts. TVie measurements are used to generate three 
tables of pressure, temperature and relative humidity at- fixed altitudes 
along the balloon trajectories- Since the altitudes at which the radiosonde 
measurements were obtained will in general be different from the fixed 
altitudes and different for each balloon release, the pressure, temperature 
and relative humidity at the fixed altitudes must be Interpolated. The 
interpolation procedure is discussed in detail in Section 3. Next the phase 
and group refractivity at f e fixed altitudes along each balloon trajectory 
are calculated. The phase and group refractivities along the ray path are 
then calculated using the interpolating formula (3-3) and equation (3-4). 

The co-latitude and longitude of the ray at the fixed altitude are determined 
by assuming the ray path is a straight line. At 10° elevation and above, 
the error introduced by this assumption is negligible. The group and phase 
refractivities at the fixed altitudes along the ray path are passed to the 
ray tracing program which calculates range (Rl^) ^ range correction and 6^ 

(see Figure 1). Tlie value of 6^ is used in Marini and Murray's formula. 

RTj^ is calculated using only the measurements obtained by the balloon which 
was released from the laser site. The Thayer method ray tracing program 
which was used by 1-arini and Murray [3] was also used to obtain our results. 
Therefore, our procedure for calculating RT^ is the same as Marini and 
Murray's. Initially, a more complex ray tracing program which required the 
0 and ((> refractivity gradients was used to calculate the range correction. 
However, the results were not substantially different from those obtained 
from the simpler and faster Thayer method program. A flowchart of the ray- 
tracing procedure is illustrated in Figure 43. 




a 


V. 



READ IN THREE RADIOSMIKS, UAVELENCTH, 
ELEVATION (E), BEGINNING (K,) AND 
FINAL (Sj) azimuth. ^ 


SURFACE 

MEASUREMEirrS 


CREATE THREE 
TABLES 


AZ - C, 


USE THREE TABLES TO CALCULATE 
REFRACTIVITY AT STANDARD 
HEIGHTS ALONG RAY PATH 


USE ONE TABLE TO CALCULATE 
REFRACTIVITY AT STANDARD 
HEIGHTS ALONG RAY PATH. 


CALL RAY TRACE 
PROGRAM 


'CALL RAY TRACE l^RT. 
I PROGRAM ^ 


AZ •= AZ +10* 


MARINI AND MURRAY'S 
FORMULA 


AZ > C, 


RT'v '* •- 1 
Rli-MM 


OUTPUT 


Figure 43. Flow chart of ^ay tracing procedure. 
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APPENDIX B 
RAY TRACE RESULTS 

The results of the ray traces through the three-dimensional and spherically 
symmetric refractivity profiles are listed on pages 67 through 144. Sheet 1 
of the data suBmar> lists the radiosonde parameters for each set of three 
balloon releases. The parameters include the release date and time and the 
altitude, azimuth angle, elevation angle and range of the balloon when the 
last measum^nt was taken. The ray trace parameters are listed on sheets 
2, 3 and 4. RT^-RT^ Is the difference between the three-dimensional and 
spherically symmetric refractivity ray traces. HM in the correction obtained 
from Marini and Hurray's formula, eq. (2-9). m-RTj^ is the difference between 
Marini and Murray's formula and the ray trace through the spherically symmetric 
profile. Azimuth selections for RT^ are from the beginning azimuth 
the final azimuth ^2 increments. Table 6 is used to calculate and 

^2 as explained in Section 5. The arc distance between release sites can be 
obtained from Table 7. 


I i d* 

i 


V. 



TABLE 6 

AZIMUTH ANGLES BETWEEN RELEASE SITES 


Release 

Release Sites | 

Sites 

51 

52 

53 

54 

55 

56 

57 

58 

51 

XX 

239° 

133° 

193° 

228° 

205° 

178° 

161° 

52 

31° 

XX 

85° 

135° 

223° 

177° 

156° 

130° 

53 

313° 

265° 

XX 

229° 

243° 

224° 

197° 

179° 

54 

13° 

315° 

49° 

XX 

259° 

215° 

172° 

114° 

55 

48° 

43° 

63° 

79° 

XX 

135° 

126° 

95° 

56 

35° 

357° 

44° 

35° 

315° 

XX 

118° 

75° 

57 

2° 

346° 

17° 

352° 

306° 

298° 

XX 

41° 

58 

314° 

310° 

359° 

294° 

275° 

255° 

229° 

XX 


TABLE 7 

ARC DISTANCE BETWEEN RELEASE SITES (KILOMETERS) 


Release 

Sites 

Release Sites 


51 

52 

53 

54 

55 

56 

57 

58 

51 

XX 

130.49 

79.97 

159.94 

242.02 

256.76 

284.95 

202.03 

52 

130.49 

XX 

181.00 

101.03 

119.95 

189.40 

248.33 

220.98 

53 

79.97 

181.00 

XX 

138.07 

242.02 

256.76 

258.85 

143.10 

54 

159.94 

101.03 

138.07 

XX 

145.22 

119.95 

151.53 

159.11 

55 

242.02 

119.95 

242.02 

145.22 

XX 

99.75 

199.94 

248.33 

56 

256.76 

189.40 

256.76 

119.95 

99.75 

XX 

105.23 

180.57 

57 

284.95 

^ 0.33 

258.85 

151.53 

199.94 

105.23 

XX 

130.49 

58 

202.03 

220.98 

143.10 

159.11 

248.33 

180.57 

13J.49 

XX 





























DATA SIDMARY 


Sheet 1: radiosonde Paraneters 

Release Date 1/30/70 



STATION 

STATION 

STATIW 

» 

S2 

55 

56 

RELEASE 

TIME 

21:30 

21:30 

21:30 

ALTITUDE 

(tan) 

15.9 

16.1 

13.1 

AZIMUTH 

76.7 

82.4 

73.0 

ELEVATION 

7.5 

7.1 

6.5 

RANGE 

(km) 

121.4 

129.4 

118.4 




.. 


RTj-RTj 


\elevation 

10 ° 

20 ° 

o 

o 

o 

o 

00 

AZIMUTH \ 





150 

- 4.23 

- 2.20 

- 1.17 

-.76 

160 

- 5.16 

- 2.47 

- 1.23 

-.76 

170 

- 6.04 

- 2.73 

- 1.30 

-.77 

180 

-6.86 

- 2.97 

- 1.36 

-.77 

190 

- 6.04 

- 2.73 

- 1.30 

-.77 

200 

- 5.16 

- 2.47 

- 1.23 

-.76 

210 

- 4.23 

- 2.20 

- 1.18 

-.76 

220 

- 3.28 

- 1.93 

- 1.10 

-.75 


SITE STN 


hW-RTj 

(cm) 


13.59 

7.04 

3.77 

2.46 

0.43 

0.39 

0.18 

0.12 






























































































































STATION 


# 


RELEASE 

TIME 


ALTITUDE 

(km) 


AZIMUTH 


ELEVATION 


52 


2:30 


15.4 


86.6 


7.8 


RANGE 

(km) 


113.6 


STATION 


STATION 




55 

58 

2:30 

2:30 

15.9 

16.0 

93.5 

84 ' 

7.0 

10.3 

130.6 

90.0 


































































RTj-RTj 




DATA SIM4ARY 




Sheet 3: 

Ray Trace Parameters 



SITE 

STN 55 

RELEASE DATE 1/27J10 

TIME 

2:30 

SITE 

ALTITUDE 

(m) 140 WAVE LENGTH 

(nm) 

530 


\ELEVATION 

10 ° 

20 ° 

o 

o 

o 

o 

00 

AZIMUTH \ 





20 

0.29 

-.10 

-.13 

-.10 

30 

0.33 

-.09 

-.12 

-.10 

40 

0.33 

1 

o 

yo 

-.12 

-.10 

50 

0.31 

1 

o 

-.12 

-.10 

60 

0.26 

-.10 

-.13 

-.10 

70 

0.18 

-.12 

-.13 

-.10 

80 

0.07 

-.14 

-.14 

-.10 

90 

-.05 

-.18 

-.14 

-.10 


SITE STN 
» 55 


MM-RT, 


13,51 
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DATA SUN4ARY 


Sheet 4 : Ray Trace Paraaeters 

SITE STN S6 RELEASE DATE 1/21/7 0 TIME 21:34 
SITE ALTITUDE («) 85 H*VE LENGTH (na) 550 



AZIMUTH \ 

10° 

o 

o 

CM 

o 

o 

80° 

RT.-Rl 
1 i 

(ca) 

340 

3.09 

m 

0.29 

m 

350 

3.17 

0.98 


0.08 

360 

3.15 

0.97 

0.29 

0.08 

370 

3.17 

0.98 

0.29 

m 

380 

3.09 

0.97 


0.08 

390 

2.94 

0.93 

0.28 

mm 

400 

2.72 

0.88 

mm 

0.08 

410 

2.42 

0.80 

0.25 

0.07 










1 

SITE ‘■'TN 
56 

m 

(m) 

13.76 

7.13 

3.81 

2.49 

(•l-RTj 

(cm) 

1.31 

0.73 

0.34 

0.23 























STATION 



ELEVATION 


STATION 


STATION 




54 

56 

25: 35 

25:51 

12.9 

13.2 

105.4 

99.5 













































































































DATA SU»*1ARY 


Sheet 4 : Ray Trace Parameters 

SITE STN 56 RELEASE DATE 1/21/7 0 TIME 23:31 
SITE ALTITUDE (m) 85 WAVE LENGTH (nm) 530 


\elevation 

10° 

o 

o 

o 

o 

o 

o 

00 

AZIMUTH \ 





340 

4.95 

1.66 

0.56 

m 

350 

4.88 

1.62 

0.55 

0.23 

360 

4.71 

1.56 

0.53 

0.22 

370 

4.88 

1.62 

0.55 

0.23 

380 

4.95 

1.66 

0.56 

0.23 

590 

4.89 

1.66 

0.56 

0.23 

400 

4.73 

1.62 

0.55 

0.23 

410 

4.46 

1.52 

0.53 

0.23 


SITE STN 
» 56 


(m) 

13.79 

7.14 

3.82 

2.50 

KM-RTj 

(cm) 

1.46 

0.64 

0.27 

0.17 

































































































































































































































STATION 



RELEASE 

TIME 

19:29 

ALTITUDE 

16.0 

(km) 

AZIMUTH 

69.8 

ELEVATION 

7.8 

RANGE 

(km) 

118.3 


STATION 


STATION 


56 

58 

19:32 

19:38 

13.5 

16.6 

70.3 

69.0 

7.8 

8.0 

100.2 

119.0 




















DATA SUMMARY 


RTj-RTj 


Sheet 2: Ray Trace Parameters 

SITE STN 52 RELEASE DATE 1/2Q/7 Q TIME i9;29 
SITE ALTITUDE fm) 85 WAVE LENGTH (m..; 530 


N^levation 

10° 

20° 

40° 

00 

0 

0 

AZIMUTH \ 





130 

-2.92 

-1.37 

-.67 

-.41 

140 

-3.27 

-1.46 

-.69 

-.41 

ISO 

-3.58 

-1.55 

-.71 

-.41 

160 

-3.85 

-1.62 

-.73 

-.41 

170 

-4.08 

-1.67 i 

-.74 

-.42 

1 

I 180 

-4.24 

1 

-1.71 

-.75 

-.42 

190 

-4.08 

-1.67 

-.74 

-.42 

200 

-3.85 

-1.62 

-.73 

-.41 

210 

-3.58 

-1.55 

-.71 

-.41 


SITE STN 


MM 

(m) 

13.60 

7.04 

3.77 

2.46 

»l-RTj 

(cm) 





0.89 

0.37 

0.14 

0.09 


























































MM-RT^ 

fan) 


1.25 




SIMIARY 

rrace Parameters 
LEASE DATE 1 / 20/7 0 TIME 15:32 
85 WAVE LENGTH (nm) 530 


20 ° 

40° 

o 

o 

00 

-.52 

-.45 

-.37 

-.49 

-.44 

-.36 

-.49 

-.44 

-.36 

-.51 

-.45 

-.36 

-.54 

-.46 

-.36 

-.59 

-.47 

-.36 

-.65 

-.48 

-.37 

-.73 

-.50 

-.37 

-.82 

-.52 

-.37 


7.0b 3.78 


0.48 


0.18 


0.11 


































































































STATION 


# 


RELEASE 

TIME 


ALTITUDE 

(km) 


AZIMUTH 


ELEVATION 


52 


17:32 


15.9 


92.5 


10.7 


STATION 


STATION 


54 

56 

'.';32 

17:30 

16.0 

16.1 


93.8 


82.2 






























































































































































































































DATA SttttARY 


9teet 3: Ray Trace Paraaeters 

SITE STN SS RELEASE DATE 1/27/7 0 TIME 6^30 






Parameters 


EASE DATE \l2ip a TIME 6:3S 
140 WAVE LENGTH (nm) 530 


o 

o 

40° 

80° 

-.32 

-.12 

-.03 

-.19 

-.09 

-.03 

-.05 

-.05 

1 

• 

O 

N) 

0.08 

1 

b 

-.02 

0.21 

0.01 

-.02 

0.32 

0.04 

1 

• 

o 

N> 

0.42 

0.06 

o 

1 

0.50 

0.08 

1 

I -.01 


0.19 


0.07 


0.04 










































f < 










































































































































































































































SlfMARY 


Trace Paraaeters 

ELEASE DATE 1/20/70 TIME 17:30 
140 WAVE LENGTH (nm) 530 


o 

o 

40° 

1 

o 

o 

00 

-.28 

-.17 

-.14 

-.17 

-.15 

-.14 

-.06 

-.12 

-.14 

0.04 

-.10 

-.14 

0.12 

-.08 

-.13 

0.19 

-.06 

-.13 

0.24 

-.05 

-.13 

0.28 

-.04 

-.13 


7.02 3.75 2.45 


0.57 


0.24 


0.15 
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DATA SUMMARY 


Sheet 2 : Ray Trace Paraneters 

SITE STN S3 RELEASE DATE 1/19/7O TIME 21:35 
SITE ALTITUDE (m) WAVE LENGTH (nn) 



\ELEVATI0N 
AZIMUTH \ 

10° 

0 

0 

<M 

0 

0 

0 

0 

00 

RTj-RTj 

(cm) 

170 

-6.04 

-2.38 

-1.02 

-.56 

180 

-6.27 

-2.45 

-1.04 

-.56 

190 

^0 

-2.38 

-1.02 

-.56 

200 

-5,72 

-2.29 


-.56 

210 

-5.31 

-2.18 

-.98 

-.56 

220 

-4.82 

-2.05 

-.94 

-.55 

230 

-4.28 

-1.90 

-.91 

-.55 

240 

-3.69 

-1.74 

-.87 

-.55 

250 

-3.08 

-1.58 

-.83 

-.54 






■ 

MM 

(m) 

13.54 

7.02 

3.75 

2.45 

KW-RTj 

(cm) 

0.84 

0.36 

0.13 

0.10 








































Sheet 3 : Ray 


SITE STN 55 REl 


Sr-E ALTITUDE (m) 



Parameters 


LEASE DATE 1/19/7 0 TIME 21:32 
140 WAVE LENGTH (run) $30 


20*^ 

40° 

00 

o 

o 

-.59 

-.59 

1 

• 

cn 

O 

-.67 

-.61 

-.50 

-.77 

-.63 

-.50 

-.90 

-.66 

-.51 

-1.04 

-.69 

-.51 

-1.20 

-.73 

-.51 

-1.37 

-.77 

-.52 

-1.53 

-.Rl 

-.52 

-1.70 

-.85 

-.52 




7.03 

3.76 

2.46 

0.32 

0,11 

0.08 






























































































































































DATA SIAMARY 


Sheet 3: Ray Trace Paraneters 

SITE STN S5 RELEASE DATE 1/19/7 0 TIME 17:32 
SITE ALTITUDE (m) 140 WAVE LEhKSTH (nm) 530 



\1LEVATI(M 

10° 

20° 

o 

o 

80° 


AZIMUTH \ 






10 

-1.71 

-1.02 

-.57 

-.38 


20 

-1.55 

-.98 

-.56 

-.38 


30 

-1.40 

-.94 

-.55 

-.37 


40 

-1.28 

-.90 

-.55 

-.37 

RTj-RTj 

(c*) 

50 

-1.18 

-.88 

-.54 

-.37 

60 

-1.10 

-.85 

-.53 

-.37 


70 

-1.06 

-.84 

-.53 

-.36 


80 

-1.05 

-.83 

-.53 

-.36 


90 

-1.07 

-.83 

-.53 

-.35 







SITE STN 

m 

(») 

13.58 

7.04 

3.76 

2.46 

# 55 

»W-RTj 

(cm) 

0.72 

0.31 

0.10 

0.07 






DATA SlMfARY 


Sheet 4 : Ray Trace Parameters 

SITE STN 58 RELEASE DATE 1/19/7 0 TIME 17:30 
SITE ALTITUDE (m) 140 WAVE LENGTH (nm) 530 


\ELEVATI0N 

10 ° 20 ° 

AZIMUTH \ 


970 -1.35 


-1.33 -.98 


-1.34 -.99 


-1.38 -1.00 


-.61 -.44 


-.62 -.43 


320 

-1.54 

-1.05 

-.63 

-.44 

330 

-1.67 

-1.08 

-.64 

-.43 

340 

-1.81 

-1.12 

-.64 

-.43 

350 

-1.97 

-1.17 

-.66 

-.44 


SITE STN 
# 58 


M4 RT, 


13.56 7.03 


0.87 0.20 


3.76 2.46 


0.16 0.01 
















































* 


RELEASE 

TIME 


ALTITUDE 


(kn) 


STATION 


52 


17; 30 


IS. 8 


AZIMUTH 


ELEVATION 


89.7 


6.6 


RANGE 

(km) 


137.5 


STATION 


STATION 






53 

56 


17;30 

17:30 


15.0 

13.1 


85.6 

82.8 


4.6 

5.9 


187.0 

127.4 
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DATA SUMMARY 


Sheet 2: Ray Trace Parameters 

SITE STN 52 RELEASE DATE 2/13/70 TIME 21;3Q 
SITE ALTIIDDE (m) 85 WAVE LENGTH (nm) 530 



\elevation 

AZDOTTH \ 

10° 

20° 

— 

40° 

o 

o 

00 

RTj-RTj 

(cm) 

100 

. 

-.55 

-.02 

0.07 

0.06 

110 

-.95 

-.13 

0.04 

0.06 

120 

-1.31 

-.23 

0.02 

wm 

130 

-1.62 

-.32 

-.00 

0.05 

140 

-1.86 

-.38 

-.02 

0.05 

150 

-2.03 

-.43 

-.03 

mm 

160 

-2.13 

-.46 

-.04 

0.05 

170 

-2.15 

-.47 

-.04 

0.04 











■ 

m 

(m) 

13.68 

7.09 

3.79 

2.48 

MM-RTj 

(cm) 

0.55 

0.35 

0.14 

0.11 


a# 
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DATA SUMMARY 


Sheet 3: Ray Trace Parameters 

SITE STN S4 RELEASE DATE 2/13/70 TIME 21:30 
SITE ALTITUDE (ra) 8S WAVE LENGTH (nm) 550 



\elevation 

AZIMUTH \ 

10° 

20° 

o 

o 

80° 

RTj.RTj 

(cm) 

230 

-1.86 

-.45 

o 

« 

1 

0.03 

240 

-1.55 

-.36 

mM 

0.03 

250 

■i 

-.27 

^^03 


260 

1 

• 

00 

o 

-.16 

o 

o 

• 

1 

0.03 

270 

-.36 

-.04 

0.02 

0.04 

280 

0.09 

0.09 

0.05 

0.04 

290 

0.54 

0.21 

0.08 

0.04 

300 

0.99 

0.33 

0.11 

0.05 











■ 

M4 

(m) 

13.65 

7.08 

3.79 

2.47 

MM-RTj 

(cm) 

0.68 

0.48 

0.22 

0.15 
















































































































































































































































































































































DATA SU»MARY 


RTj-RTj 


Sheet 3: Ray Trace Paraneters 

SITE STN S3 RELEASE DATE 1/26/70 TIME 9:30 

SITE ALTITUDE (m) 85 WAVE LENGTH (nn) 530 


\elevation 
AZIMUTH X 


20° 

o 

o 

o 

o 

00 

-1.77 

-.73 

-.37 

-1.64 

1 

o 

• 

1 

-.37 

-1.49 

-.67 

-.36 

-1.34 

-.63 

-.36 

-1.17 

-.59 

-.36 

-.99 

-.55 

-.35 

-.82 

-.51 

-.35 

-.66 

-.47 

-.34 


SITE STN 


13.42 6.96 5.73 


NW-RT^ 

(cm) 


-.11 


-.02 


-.06 


-.04 





























































































































































STATION 


* 


RELEASE 

TIME 


ALTITUDE 

(km) 


AZIMUTH 


ELEVATION 


52 


21:32 


15.9 


106.1 


7.2 


RANGE 

(km) 


127.1 


SUIMARY 


T 


iiosonde Parameters 
ate 1/21/70 



STATION 

STATION 


55 

56 


21:47 

21:34 


13.2 

12.6 


104.3 

96.9 


6.4 

5.1 





119.1 


142.1 

































































DATA SIM4ARY 


Sheet 1: Radiosonde Paraneters 

Release Date 1/21/70 



STATION 

STATION 

STATION 

# 

52 

54 

56 

RELEASE 

TIME 

23:36 

23:35 

23:31 

ALTITUDE 

(km) 

15.9 

12.9 

13.2 

AZIMUTH 

107.4 

105.4 

99.3 

ELEVATION 

7.7 

4.9 

6.3 

RANGE 

(km) 

118.4 

149.4 

120.4 


DATA SUKMARY 


RTj-RTj 


Sheet 2: Ray Trace Parameters 

SITE STN 56 RELEASE DATE 1/21/70 TIME 23:31 
SITE ALTITUDE (m) 85 WAVE LENGTH (nm) 353 


\ELEVATI0N 

A2IMUTHX 


10° 

s 

20° 

o 

o 

o 

o 

00 

5.35 

1.79 

0.61 

m 

5.29 

1.76 

0.60 

0.24 

5.10 

1.68 

0.58 

0.24 

5.29 

1.76 

0.60 

0.23 

5.35 

1.79 

0.61 

tm 

5.30 

1.79 

0.61 

0.23 

5.12 

1.76 

0.60 

0.23 

4.83 

1.69 

0.59 

0.23 


SITE STN 
# 56 


M4 

(m) 

14.91 

7.73 

4.13 

2.70 

MM-RT^ 

(cm) 

1.46 

0.59 

0.25 

0.17 
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DATA SUMMARY 


Sheet 2: Ray Trace Parameters 

SITE STN 56 RELEASE DATE 1/30/7 0 TIME 21:30 
SITE ALTITUDE (m) 140 WAVE LENGTH (nm) 353 



\EL£VATION 
AZIMUTH \ 

10° 

20° 

o 

o 

o 

o 

00 

RTj-RTj 

(cm) 

300 

0.97 

-.97 

-1.05 

-.93 

310 

0.97 

-.98 

-1.05 

-.92 

320 

0.79 

-1.03 

-1.06 

-.93 

330 

0.43 

-1.14 

-1.09 

-.93 

340 

o 

1 

-1.29 

-1.12 

-.93 

350 

1 

-.75 

-1.49 

-1.17 

-.94 

360 

-1.54 

-1.72 

-1.23 

-.94 

370 

-.75 

-1.49 

-1.17 

-.94 











SITE STN 
# 56 

MM 

(m) 

14.67 

7.60 

4.07 

2.66 

(cm) 

0.75 

0.52 

0.25 

0.18 
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